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ABSTRACT 

We present mid-infrared imaging at 11.7 fim for the quadruple lens systems, PG1115-I-080 and 
B1422-I-231, using the cooled mid-infrared camera and spectrometer (COMICS) attached on the Sub¬ 
aru telescope. These lensed QSOs are characterized by their anomalous optical and radio flux ratios as 
obtained for (Al, A2) images of PG1115-I-080 and (A, B, C) images of B1422-I-231, respectively, i.e., 
such flux ratios are hardly reproduced by lens models with smooth mass distribution. Our mid-infrared 
observations for these images have revealed that the mid-infrared flux ratio A2/A1 of PG1115-I-080 
is virtually consistent with smooth lens models (but inconsistent with the optical flux ratio), whereas 
for B1422-I-231, the mid-infrared flux ratios among (A, B, C) are in good agreement with the radio 
flux ratios. We also identify a clear infrared bump in the spectral energy distributions of these QSOs, 
thereby indicating that the observed mid-infrared fluxes originate from a hot dust torus around a QSO 
nucleus. Based on the size estimate of the dust torus, we place limits on the mass of a substructure in 
these lens systems, causing the anomalous optical or radio flux ratios. For PG1115-I-080, the mass of 
a substructure inside an Einstein radius, Me, is < 16Mq, corresponding to either a star or a low-mass 
CDM subhalo having the mass of Mf^Q < 2.2 x inside radius of 100 pc if modeled as a singular 

isothermal sphere (SIS). For B1422-I-231, we obtain Me ^ 2O9M0, indicating that a CDM subhalo is 
more likely, having the mass of Mf^Q > 7.4 x IO^'^Mq. 

Subject headings: gravitational leasing — infrared: galaxies — quasars: individual (PG1115-I-080, 
B1422-f231) 


1. INTRODUCTION 

The cold dark matter (CDM) scenario for structure 
formation in the Universe has been quite successful to 
explain a wide variety of observational results, such as 
the characteristic fluctuations of the cosmic microwave 
background (CMB) and large-scale structure of galaxy 
distribution, on the spatial scales larger than ~ 1 Mpc. 
In particular, the recently released WMAP results in¬ 
dicate that the power spectrum of the CMB radiation 
matches the CDM prediction remarkably well (Spergel 
et al. 2003). The currently best world model, which ac¬ 
cords with the results of many cosmological observations, 
consists of approximately 23 % CDM, 4 % baryon, and 
73 % vacuum energy. 

However, the recent advent of high-resolution N-body 
simulations on CDM-based structure formation has en¬ 
abled to highlight various discrepancies with existing ob¬ 
servations on the spatial scales smaller than ^ 1 Mpc. 
One of the most serious issues is that CDM models pre¬ 
dict the existence of several hundred dark satellites or 
“CDM subhalos” (with masses of M ^ 10^“® Mq) in a 
galaxy-sized halo (with M ~ 10^^ Mq), in sharp con- 

^ Based on data collected at Subaru Telescope, which is oper- 
ated by the National Astronomical Observatory of Japan. 

^ Astronomical Institute, Tohoku University, Aoba-ku, Sendai 
980-8578, Japan; chiba@astr.tohoku.ac.jp 

^ Institute of Astronomy, School of Science, University of Tokyo, 
Mitaka, Tokyo 181-0015, Japan; minezaki@mtk.ioa.s.u-tokyo.ac.jp 

^ National Astronomical Observatory, Mitaka, Tokyo 181-8588, 
Japan; kashik@optik.mtk.nao.ac.jp 

® Institute of Space and Astronautical Science, Japan Aerospace 
Exploration Agency, Sagamihara, Kanagawa 229-8510, Japan; 
kataza@ir.isas.jaxa.jp 

^ School of Science and Engineering, Kinki University, Higashi 
Osaka 577-8502, Japan; kinoue@phys.kindai.ac.jp 


trast to the observed number of about a dozen Milky 
Way satellites (e.g., Klypin et al. 1999; Moore et al. 

1999) . Although this discrepancy may be alleviated by 
some baryonic processes, such as the suppression effect 
of background UV radiation on the formation of visi¬ 
ble dwarf galaxies (e.g., Bullock, Kravtsov, & Weinberg 

2000) , the presence of many dark subhalos in a galaxy is 
inevitable if CDM models are correct. It is worth not¬ 
ing that the abundance of CDM subhalos is closely rele¬ 
vant to the small-scale power spectrum of initial density 
fluctuations, which determines the formation of the first 
(Population HI) stars and resultant reionization history 
of the Universe (e.g., Yoshida et al. 2003). Also, the 
abundance of such subhalos may be related to the for¬ 
mation of bright galaxies like the Milky Way, as their 
clustering history seems to affect the final morphologi¬ 
cal type of forming galaxies (e.g., Kauffmann, White, & 
Guiderdoni 1993). Thus, clarifying the presence or ab¬ 
sence of CDM subhalos in a galaxy like the Milky Way is 
closely related to how visible parts of the Universe have 
formed. 

Metcalf & Madau (2001) and Chiba (2002) demon¬ 
strated that gravitational leasing is the most power¬ 
ful probe of such numerous CDM subhalos which re¬ 
side in external lensing galaxies. They showed that if 
these subhalos exist within lensing galaxies responsible 
for multiple imaged QSOs, the magnification properties 
of lensed images are modified compared to the predic¬ 
tions of smooth lens models - as was first noticed by 
Mao & Schneider (1998) who examined the effects of 
non-smooth density distribution in lenses. In particu¬ 
lar, the anomalous image flux ratios observed in some 
four-image QSOs, which are hardly reproduced by any 
lens models with smooth density distribution, can be ex- 
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plained by the existence of numerous CDM subhalos in 
the lensing galaxies, while the image positions are little 
affected (Chiba 2002). More extensive, statistical anal¬ 
yses of anomalous image flux ratios have provided im¬ 
portant limits on the nature of lens substructures and 
also the power spectrum of CDM models at small spa¬ 
tial scales (e.g., Metcalf & Zhao 2002; Dalai & Kochanek 
2002; Keeton 2003; Kochanek & Dalai 2004; Oguri 2004). 
However, such a lensing approach to probing numerous 
CDM subhalos is subject to generic uncertainties asso¬ 
ciated with the microlensing by stellar populations in 
lensing galaxies, which causes the same magnification ef¬ 
fect on the QSO images (Schechter & Wambsganss 2002). 
Thus, to map out the CDM subhalos alone, we require 
more convincing observational studies in addition to the 
existing image data. 

One of the best techniques to distinguish microlens¬ 
ing by stars from lensing by CDM subhalos"^ is to ob¬ 
serve flux ratios from emission regions of different size, 
Rs, and to compare them with the mass-dependent 
size of the Einstein radius, which is, for a point-mass 
lens model, given as Re{M) ~ O.Ol(M/M 0 )^/^/i“^/^ pc 
(where h = Hq/IOO km s“^ Mpc“^) at the typical red- 
shift of a QSO at zs = 2—3 and a lens at zl = 0.3—0.5. If 
Rs ^ Re{M) for either a star or subhalo with mass M, 
then the image is unaffected by such a lens perturber, 
whereas it is significantly magnified if Rs < Re{M). 
This indicates that the observed optical continuum emit¬ 
ted from a QSO nucleus, a central engine of an accretion 
disk with an emission region of Rs ~ 10 ^^ cm (e.g., 
Wambsganss, Schneider, & Paczyhski 1990; Rauch & 
Blandford 1991; Wyithe et al. 2000), is subject to lens¬ 
ing magnification by both stars and CDM subhalos. On 
the other hand, if we focus on spatially extended regions 
around the central engine, such as a surrounding dust 
torus responsible for infrared emission of a QSO (e.g., 
Sanders et al. 1989) or line-emitting regions (e.g., Mous- 
takas & Metcalf 2003), then stars would be unable to 
affect the flux from such extended regions of i?s > 1 pc 
(Agol, Jones, & Blaes 2000; Metcalf et al. 2004). Thus, 
this selective magnification offers an important probe for 
distinguishing the effects of CDM subhalos from those of 
stars on the concerned QSO images showing anomalous 
optical flux ratios. 

In this paper, we report the mid-infrared observations 
of lensed QSOs showing anomalous optical flux ratios 
using the Subaru telescope. The observed mid-infrared 
waveband corresponds to the near-infrared waveband in 
the rest frame, and its flux is considered to be dominated 
by thermal radiation from hot dust located at the inner¬ 
most region of a dust torus (Barvainis 1987; Kobayashi 
et al. 1993). The inner radius of a dust torus would 
be determined by the highest sublimation temperature 
of dust (T ~ 1800 K) and the UV luminosity of a QSO 
central engine, and it is estimated as Rs ~ 1 pc for 
UV luminosity of 10'^® erg s“^ (Barvainis 1987). There¬ 
fore, the observed mid-infrared flux is totally insensitive 
to microlensing by foreground stars {Rs ^ Re) but is 
strongly magnified by CDM subhalos {Rs Re)- 

Moreover, mid-infrared observations of lensed images 
are advantageous compared to other wavebands in this 

^ This lensing is often referred to as millilensing because its 
Einstein radius is generally of the order of milliarcsec. 


kind of study: ( 1 ) mid-infrared flux is almost free from 
the effects of differential extinction among different im¬ 
ages, in contrast to the use of continuum or line fluxes in 
shorter wavebands, ( 2 ) imaging data of both radio-quiet 
and radio-loud QSOs are available, and (3) the source 
size of mid-infrared flux, which originates from near- 
infrared emission in the rest frame, is quantitatively es¬ 
timated from the rest-frame optical luminosity of a QSO 
based on the result of the dust reverberation observa¬ 
tion, i.e. the measurement of arrival time between UV 
flux from a central engine and near-infrared flux from a 
surrounding dust torus (Minezaki et al. 2004). The last 
point is of particular importance to place useful limits 
on the mass of substructure in a foreground gravitational 
lens. 

The paper is organized as follows. In §2, we show the 
target selection and the observations. In §3, the proce¬ 
dure of data reduction is presented to obtain the mid- 
infrared fluxes and their ratios among the lensed images. 
In §4, the implications for the current observational re¬ 
sults are discussed and the conclusions are drawn in §5. 
In what follows, we adopt the set of cosmological param¬ 
eters of H = 0.3, A = 0.7, and h = 0.7 for all relevant 
estimations. 

2. OBSERVATION 
2 . 1 . Targets 

We select two lens systems with four images, 
PG1115-b080 at zs = 1.72 and B1422-b231 at zs = 3.62, 
renowned for their anomalous image flux ratios. The 
lensing galaxies are possibly elliptical galaxies at Zl = 
0.31 and 0.34, respectively (Kundic et al. 1997a; Kundic 
et al. 1997b; Tonry 1998). The former lens system holds 
the closely separated pair of images Al and A2 with a 
separation of 0."48, and this configuration emerges if 
the QSO is close to and inside a fold caustic provided 
by the foreground lens. The latter shows the colinear, 
three highly magnified images. A, B, and C, and this 
configuration emerges if the QSO is close to and inside 
a cusp caustic. In such lens systems associated with a 
fold or cusp caustic, there exists a universal relation be¬ 
tween the image fluxes, i.e., A2/A1= 1 or (A-|-C)/B= 1, 
whereas the observed flux ratios violate these rules sig¬ 
nificantly, A2/A1= 0.64 ± 0.02 (Impey et al. 1998) and 
(A-|-C)/B= 1.50 ±0.01 (Patnaik et al. 1992). It is worth 
noting that the flux ratio of A2/A1 in PG1115±080 
showed little variation with wavelength from the multi¬ 
ple wavelength observations by Impey et al. (1998), and 
that the flux ratios observed in B1422±231 by Patnaik 
et al. (1992) are based on radio fluxes, so the observed 
flux ratios of the concerned images are unlikely to be af¬ 
fected by dust extinction. Thus, these two lenses are the 
most representative targets for investigating the effects 
of substructure in lensing galaxies. 

2 . 2 . Observation 

The mid-infrared imaging of these lensed QSOs was 
carried out on the nights of UT 2004 May 5 and 6 , 
using the cooled mid-infrared camera and spectrome¬ 
ter (GOMICS; Kataza et al. 2000) attached on the 
Gassegrain focus of the Subaru telescope. The field of 
view is 42" x 32" and the pixel scale is 0."129 pixel“^. 
We used the Nil.7 filter, whose effective wavelength and 
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Fig. 1.— The quadruple lens systems, PG1115+080 (left) and B1422+231 (right), at 11.7 pm taken with COMICS/Subaru on 2004 
May 5 and 6 (UT). The direction of the images is that the north is up and the east is left, and the pixel scale of them is 0.”065 pixel“^. 
These images have been smoothed with a Gaussian kernel of tr = 0.”065 in order to improve their visual impression. 


bandwidth are Ac = 11.67 fj,m and AA = 1.05 /im, re¬ 
spectively. The chopping frequency was set 0.45 Hz with 
a width of 10", then the target images were within the 
field of view for both chopped positions. In addition 
to the chop, the telescope position was nodded several 
times during the observation. It was photometric except 
for the first quarter of both nights, and the diffraction 
core of point-spread function (PSF) could be seen. The 
FWHM of PSF was 0."33 at small airmass. Since we 
found some degradation of PSF for the images at large 
airmass, we limited the data to those at small airmass, 
< 1.6 for PG1115-b080 and < 1.4 for B1422-b231. The 
total exposure times of the available data were 1.8 hours 
for PG1115-I-080 and 3.1 hours for B1422-I-231, respec¬ 
tively. We also observed HD98118 and HD127093 for the 
photometric standard stars (Gohen et al. 1999). 

The images were reduced using IRAF ®. An image data 
file consisted of 208 successive chopped frames and its ex¬ 
posure time was 200.5 seconds. First, the sky background 
of the frame of one chopped position was subtracted by 
the frame of the other chopped position, and flat-field 
correction was applied using a sky flat image that was 
assembled by combining the frames of the other chopped 
position and correcting global variation across the im¬ 
age. Then two images with exposure time of 100 seconds 
were obtained from an image data file by combining the 
reduced frames of each chopped position. After that, in 
order to reduce residual sky background around the tar¬ 
gets, we fitted it with a low order function on the square 
area of 8.3" x 8.3" around the target (masking center 

® IRAF is distributed by the National Optical Astronomy Ob¬ 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 


3.1" X 3.1") and subtracted it. The small systematic drift 
of image positions due to the telescope tracking error was 
fitted with smooth functions using the image positions of 
QSO themselves which were detected marginally in each 
100 seconds-exposure image, then those images were reg¬ 
istered according to the fitted drift functions. In order 
not to degrade the angular resolution by the sub-pixel 
shifts during the image registration, drizzle (Fruchter & 
Hook 2002) task was used to perform sub-pixel shifts and 
re-sample with smaller size of pixels than the original, 
where the task parameters were set as p = 1.0 (equiva¬ 
lent to shift-and-add) and s = 0.5 (re-sampled with half 
size of pixels). Finally, all available data were combined. 

The resultant mid-infrared (A = 11.7 pm) images 
of the lensed QSOs, PG1115-I-080 and B1422-I-231, are 
presented in Figure 1. The images presented in the 
figure have been smoothed with a Gaussian kernel of 
cr = 0."065 in order to improve their visual impression, 
although photometries were carried out based on the im¬ 
ages without any smoothing. As presented in Figure 1, 
the lensed images Al and A2 of PG1115-I-080 and A, B, 
and C of B1422-1-231 were clearly detected and well sepa¬ 
rated from each other. The minimum image separations, 
0."48 and 0."50 for these lens systems, respectively, are 
3.0 and 3.1 times larger than ar of PSF (cr^. = 0."16, as 
will be described later). The faint lensed images B and 
G of PGl 115-1-080 were also detected. 

3. RESULTS 
3.1. Flux Ratio 

The flux ratios between the lensed images of these 
QSOs were estimated by PSF fitting photometry. We 
assumed a circular, Gaussian radial profile for the PSF 
model, and the relative positions between the images 
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TABLE 1 

Mid-infrared Flux Ratio and Flux of PG1115-I-080 


TABLE 3 

Comparison of Flux Ratios for PG1115+080 



flux ratio 


flux 


A2/A1 

B/Al 

C/Al 

A1±A2 

(mjy) 

All 

(mJy) 

0.93 ±0.06 

0.16 ± 0.07 

0.21 ±0.04 

14.6 ±1.2 

17.5" 


“^The total flux of all four images, based on the flux of A1-I-A2 
and the flux ratios listed here. 


TABLE 2 

Mid-infrared Flux Ratio and Flux of B1422+231 



flux ratio 


flux 

(A±C)/B 

A/B 

C/B 

A±B±C 

(mJy) 

1.51 ±0.06 

0.94 ±0.05 

0.57 ±0.06 

19.2 ± 2.9 


were taken from the CASTLES web site®, which are 
based on imaging observation using the Hubble Space 
Telescope (HST). The free parameters for this fit were 
the fluxes of the lensed images, the positional shift of the 
whole images, and Ur for the Gaussian radial profile of 
PSF. We used (j) = l."3 around each lensed image for the 
fitting aperture, although the estimated flux ratios were 
almost independent of the size of the fitting aperture; 
using a Gaussian radial profile for the model PSF, the 
diffraction core of the real PSF was fitted effectively. In 
the case of PG1115+080, we first fit only the images A1 
and A2, then with the derived parameters being fixed we 
fit the images Al, A2, B, and C (where the fluxes of B 
and G are free parameters). In the case of BI422-I-231, 
we fit only the images A, B, and G because the faint 
image D was not detected. 

Since the error of flux ratio was not determined only 
by photon statistics and detector noise, and since the 
fluctuation of sky background was important, the error 
of flux ratio was estimated by the following simulation. 
First, eight blank-sky areas of 8.3" x 8.3" were selected 
surrounding the objects, and the model images of targets 
were added to the sky areas of the reduced image with ex¬ 
posure time of 100 seconds. The simulated model image 
has a Gaussian radial profile with the best fit parameters 
of the total flux and CFr- The flux ratio of the simulated 
model images was also the same as that of the best-fit to 
the observation. Then, those model images were reduced 
in the same way as on the real images, and the PSF fit¬ 
ting photometries were applied to them to measure flux 
ratios. Finally, the error of the observed flux ratio was 
derived from the standard deviation of the flux ratios of 
the eight simulated images. 

The resultant mid-infrared flux ratios and their one- 
sigma errors are presented in Table 1 and 2. The best- 
fit values of ar to the observation were 0."16 for both 
lensed QSOs. For comparison. Table 3 and 4 show the 
prediction of a smooth lens model (Ghiba 2002) and the 
flux ratios in other wavebands reported for the concerned 
images^®. It follows that (1) the mid-infrared flux ra- 

® http://cfa-www.harvard.edu/glensdata/ 

Errors listed in these tables are, when available, based on the 


Wavelength 

Date 

A2/A1 

B/Al 

C/Al 

Lens Model" 


0.92 

0.22 

0.28 

11.7 ^m'’ 

2004.05.05-06 

0.93 ± 0.06 

0.16 ±0.07 

0.21 ± 0.04 

F160W" 

1997.11.17 

0.63 ± 0.02 

0.16 ±0.01 

0.25 ± 0.01 

F814W" 

1997.05.17 

0.67 ± 0.02 

0.17 ±0.01 

0.25 ± 0.01 

F555W" 

1999.03.31 

0.52 ±0.07 

0.15 ±0.02 

0.25 ± 0.04 

F160W“' 

1997.11.17 

0.64 ± 0.02 

0.17 ±0.01 

0.26 ± 0.01 

F160W® 

1997.11.17 

0.67 ± 0.04 

0.22 ±0.02 

0.29 ± 0.01 

K’f 

1999.01.11-13 

0.67 

0.19 

0.28 

Jf 

1999.01.11-13 

0.59 

0.16 

0.24 

F785LPS 

1993.02.18 

0.71 

0.18 

0.27 

F555WS 

1993.02.18 

0.63 

0.17 

0.23 

F785LP'' 

1991.03.03 

0.70 ±0.01 

0.16 ±0.003 

0.26 ± 0.01 

F555W'’ 

1991.03.03 

0.66 ± 0.03 

0.16 ±0.01 

0.26 ± 0.01 


“‘Prediction of the smooth lens model by Chiba (2002). 

'“’This work. 

““Taken from the CASTLES web site. The observing date is taken 
from the HST data archive. 

“'impey et al. (1998). 

®Treu & Koopmans (2002). 

^Iwamuro et al. (2000). 

^Unpublished WFPCl images analyzed by Impey et al. (1998). 
''Kristian et al. (1993). 


TABLE 4 

Comparison of Flux Ratios for B 1422-1-231 


Wavelength 

Date 

(A±C)/B 

A/B 

C/B 

Lens Model" 


1.25 

0.76 

0.50 

11.7 Atm'’ 

2004.05.05-06 

1.51 ±0.06 

0.94 ±0.05 

0.57 ± 0.06 

F160W" 

1998.02.27 

1.43 ±0.05 

0.90 ±0.03 

0.53 ± 0.02 

F791W" 

1999.02.06 

1.57 ±0.13 

0.97 ±0.10 

0.60 ± 0.04 

F555W" 

1999.02.06 

1.57 ±0.18 

1.02 ±0.14 

0.55 ± 0.06 

F480LP“' 

1995.04.20 

1.34 ±0.04 

0.79 ± 0.02 

0.54 ± 0.02 

F342W“* 

1995.04.21 

1.21 ±0.04 

0.68 ±0.02 

0.53 ± 0.02 

Gunn r^ 

1993.04.26 

1.23 ±0.02 

0.75 ±0.015 

0.48 ± 0.010 

Gunn g® 

1993.04.26 

1.25 ±0.02 

0.77 ±0.020 

0.48 ± 0.010 

8.4 GHz' 

1997.06.11-12 

1.42 ±0.02 

0.93 ±0.02 

0.49 ± 0.01 

5 GHzS 

1991.08.31 

1.50 ±0.01 

0.98 ±0.01 

0.52 ± 0.01 

8.4 GHzS 

1991.06.16 

1.48 ±0.01 

0.97 ±0.01 

0.52 ± 0.01 


"Prediction of the smooth lens model by Chiba (2002). 

'’This work. 

"Taken from the CASTLES web site. The observing date is taken 
from the HST data archive. 

“'impey et al. (1996). 

““Yee & Ellingson (1994). 

'Patnaik et al. (1999). 
sPatnaik et al. (1992). 


tio, riR, obtained for A2/A1 of PG1115-I-080 is consis¬ 
tent with the prediction of a smooth lens model, rsmooth 
(~ 1), whereas the optical flux ratio, Vopt (— 0.65), is 
systematically smaller than rsmooth- (2) the mid-infrared 
flux ratios for (A, B, C) images of B1422-I-231 are virtu¬ 
ally consistent with the reported radio flux ratios (Pat- 
naik et al. 1992), i.e., the mid-infrared flux ratios re¬ 
main anomalous as well compared with the prediction of 
a smooth lens model. These properties in the flux ratios 
of both lens systems provide important implications for 
the nature of lens substructure as discussed later. 

reported magnitude errors in lensed images or relative errors in 
flux ratios. 
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3.2. Flux and Spectral Energy Distribution 

The fluxes of these lensed QSOs were estimated by 
aperture photometry. We used (f) = l."3 for the photo¬ 
metric aperture. Since the lensed images were so close 
each other that their flux distributions were overlapped 
each other, we measured the fluxes of /(A1-I-A2) of 
PG1115-f080 and f{A+B+C) of B1422-h231 by connect¬ 
ing the photometric apertures. The B and C images of 
PG1115-I-080 were so faint that aperture photometries 
for them were not executed. Instead, the total flux of all 
images of PGI115-I-080 was estimated from /(AI-I-A2) 
and the flux ratios we obtained. After measuring the 
fluxes within the apertures, the aperture corrections were 
applied to estimate the total fluxes. The amount of cor¬ 
rections was estimated as about 10%. The error of flux 
was estimated in the same way as on the flux ratio. 

The fluxes of the lensed QSOs were calibrated by 
comparing the fluxes of standard stars, HD98118 and 
HD127093. Since HD127093 is a variable star of irreg¬ 
ular type, it was used only for the estimation of an air- 
mass gradient. The flux of HD98118 was calculated by 
integrating its spectral energy distribution (SED) pre¬ 
sented by Cohen et al. (1999) with the top-hat function 
of Ac = 11.67 /im and AA = 1.05 /xm. 

The resultant fluxes and their one-sigma errors are also 
listed in Table 1 and 2, and our results of mid-infrared 
fluxes of PG1115-I-080 and B1422-I-231 are graphically 
presented in Figure 2, where fluxes in other wavebands 
are also shown for comparison (CASTLES; Barvainis & 
Ivison 2002; Iwamuro et al. 2000; Impey et al. 1998; 
Courbin et al. 1997; Christian, Crabtree & Waddell 
1987; Impey et al. 1996; Remy et al. 1993; Yee & Bech- 
told 1996; Yee & Ellingson 1994; Oyabu et al. 2001; Pat- 
naik & Narasimha 2001; Patnaik et al. 1992). It follows 
that our mid-infrared fluxes of both QSOs present clear 
infrared bumps in their SEDs, thereby indicating that 
the observed mid-infrared fluxes of PG1115-I-080 and 
B1422-I-231 are dominated by thermal emission of a dust 
torus surrounding the central black hole and accretion 
disk. The SEDs in Figure 2 are compared with the mean 
energy distribution (MED) of low-redshift QSOs by Elvis 
et al. (1994), and our mid-infrared fluxes are consistent 
with the MED, if the dispersion of SED around the MED 
is taken into account. Oyabu et al. (2001) measured the 
total flux of B1422-I-231 in mid-infrared bands using ISO. 
The observation date was UT 1997 January 13, and the 
reported fluxes were 5.8 mJy in LW2(A = 6.7 /xm) and 
15.1 mJy in LW3(A = 14.3 /xm). These fluxes seem to be 
somewhat fainter than our result, however these obser¬ 
vations are not inconsistent if the error of our total flux 
and the absolute photometric calibration error of ISO 
observation (15% in Oyabu et al. 2001) are considered. 
Therefore, mid-infrared flux variation was not found be¬ 
tween two observations of B1422-I-231. 

4. DISCUSSION 

4.1. Limits on Substructure Sensing 

The mid-infrared flux ratios presented here set impor¬ 
tant limits on the nature of lens substructure, because 
the observed flux is dominated by thermal emission of 
hot dust at the innermost region of a dust torus, which 
is extended and whose size can be estimated based on the 
results of the reverberation observations. Using the lag 
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Fig. 2.— The spectral energy distributions of PG1115-I-080 (left) 
and B1422-I-231 (right). The open circles are taken from literature, 
and the filledcircles are the 11.7 /im fluxes of this work. The dashed 
line is the mean energy distribution (MED; Elvis et al. 1994) of 
radio quiet QSOs, and the dash-dotted line is the MED of radio 
loud QSOs, both of which are normalized by the F160W fluxes 
taken from the CASTLES. 




time At from the UV-optical continuum flux variation to 
the near-infrared flux variation, which is regarded as the 
light-travel time from the central accretion disk to the 
innermost dust torus, the radius of a near-infrared emis¬ 
sion region is estimated as R = cAt where c is the speed 
of light. More luminous active galactic nuclei (AGN) 
have larger dust tori, and the lag time of near-infrared 
(iL-band, A = 2.2 /xm) flux variation can be estimated 
as log At(days) = —2.15 — My/5.0 where My is the ab¬ 
solute magnitude of AGN in V band, according to the 
results of dust reverberation observations compiled by 
Minezaki et al. (2004). 

Based on the total apparent magnitudes of the quadru¬ 
ple images in F160W taken from the CASTLES web site 
(14.9 mag for PG1115-b080 and 13.3 mag for B1422-k231) 
and the magnification factors due to lensing galaxies 
(32.73 and 22.79, respectively, estimated from the lens 
model of Chiba 2002), we obtain My = —24.6 mag and 
—28.3 mag, respectively. Then, the radii of the rest- 
frame AT-band emission regions are estimated as R = 0.5 
pc and 2.7 pc for PG1115-b080 and B1422-k231, re¬ 
spectively. The rest-frame wavelength of the observing 
11.7 ^m of PG1115-I-080 is 4.3 ^m, which is somewhat 
longer than that of K band. Concerning the contribu- 





6 


CHIBA ET AL. 


tion to the observed mid-infrared flux of cooler dust, we 
estimate the source radius as Rs ~ 1 pc. The source ra¬ 
dius of observed mid-infrared flux of B1422-I-231 is esti¬ 
mated as Rs — 2.7 pc, because its observing wavelength 
of 11.7 /im corresponds to the rest-frame wavelength of 
2.5 /im, almost the same as that of K band. At last, we 
convert the linear dimension Rs of a source image into 
an angular size of 9s, and we obtain dg ~ 1 x 10“^ arcsec 
and 3.7 x IQ-* arcsec for PG1115-b080 and B1422-f231, 
respectively. 

Next, the characteristic sizes of UV-optical continuum 
emission regions are estimated. They would be depen¬ 
dent on QSO luminosity, and also dependent on ob¬ 
serving wavelength as larger source size is expected at 
longer wavelength. We estimate the source sizes for the 
F160W observations of PG1115-I-080 and B1422-I-231 as 
i?s ~ 6 X 10“^ pc and 2 x 10“^ pc, yielding < 7 x 10“^ 
arcsec and 3 x 10“® arcsec, respectively (see Appendix 
for details). 

Either a star with sub-solar mass (microlensing) or a 
GDM subhalo with M = (millilensing) is re¬ 

sponsible for the reported anomalous flux ratios. To dis¬ 
tinguish between these substructure lensing effects, it is 
worth estimating an Einstein angle to compare with 
6s- Denoting Me as a. substructure mass inside an Ein¬ 
stein angle, the latter is expressed as, 

9e = 8x 10-^ (o;^) arcsec , (1) 

for both lens systems. This equation applies to both a 
point-mass lens and spatially extended one, provided a 
lens is circularly symmetric in projection. As a reference 
for the total mass of a spatially extended substructure, 
being applicable to a GDM subhalo, we assume a singular 
isothermal sphere (SIS), which is parameterized by a one¬ 
dimensional velocity dispersion, (j||, or circular velocity, 

Vc = ^0\\, so that the Einstein angle is written as 9 e — 
lx 10“^(14/km s~^)^ arcsec for both lens systems. Then, 
the total mass of an SIS substructure inside radius of 
100 pc is given as = 2.3 x 10"‘(14/km s“^)^ Mq. 

Compared with 9s for a surrounding dust torus or 
continuum-emitting region, one can immediately draw 
the conclusion that a star with sub-solar mass is unable 
to magnify mid-infrared flux arisen from a dust torus. 
This is also the case for low-mass GDM subhalos, pro¬ 
vided 9e «C 9s- Then, for PG1115-I-080, the reported 
anomalous flux ratios in the optical band may be caused 
by either microlensing of a star with sub-solar mass or by 
a very low-mass GDM subhalo with Me ^ 1.6 x IO^Mq 
from the 9e ^ 9s ansatz, where both of these sub¬ 
structure lenses do not magnify mid-infrared flux from 
a dust torus. Tighter mass limits are available, based 
on the detailed simulations by Wyithe, Agol, & Fluke 
(2002), which showed that 9e must be at least an or¬ 
der of magnitude smaller than 9s for no magnification. 
We then obtain Me IQMq inside an Einstein radius 
and M^qq < 2.2 x lO'^M© inside radius of 100 pc if a 
substructure is modeled as an SIS. For B1422-I-231, both 
of the radio and mid-infrared flux ratios appear to be 
modified by substructure lensing compared with the pre¬ 
diction of a smooth lens model; the difference in the most 
problematic flux ratio, A/B, is about 20 %. Taking into 
account 9s of a dust torus, the simple ansatz 9e ^ 9s 


for the amplification of a torus suggests a substructure 
with Me ^ 2.1 x IO'^Mq. Recently, Inoue & Chiba 
(2004) investigated the detailed properties of extended 
source effects in substructure lensing and showed that 
even if a source size is about ten times larger than an 
Einstein radius of an SIS-form substructure, the flux ra¬ 
tio is changed by about 20 % from a smooth-lens predic¬ 
tion, provided a substructure is just centered at a source 
image. Thus, we obtain > 209M© inside an Ein¬ 
stein radius, whereby a GDM subhalo is more likely, and 
Miqq > 7.4 X lO'^M© inside radius of 100 pc if a sub¬ 
structure is modeled as an SIS. 

4.2. Effects of Other Possible Mid-Infrared Sources 

Although the mid-infrared fluxes of both target QSOs 
appear to be dominated by thermal emission of a sur¬ 
rounding dust torus, minor fraction of them might be 
contributed by the flux of more compact regions than a 
dust torus, such as an extension to infrared wavebands of 
optical continuum emission from a central accretion disk 
(Collin & Bure 2001) or non-thermal emission from a 
compact region relating to radio activity (Neugebauer & 
Matthews 1999; Enya et al. 2002). Also, more spatially 
extended sources than a dust torus, such as a region of 
nuclear starburst and/or QSO’s host galaxy as a whole, 
are to be taken into account. We discuss here the effects 
of possible contribution to mid-infrared fluxes by both 
more compact and extended regions than a dust torus. 

4.2.1. Effects of More Compact Sources 

The contributions from the accretion disk to the mid- 
infrared flux are estimated as ~ 0.06 for PG1115-I-080 
and ~ 0.2 for B1422-I-231, by extrapolating the optical 
power-law continuum from the CASTLES F160W fluxes 
with index a,^ = 0 where oc 1 /“" [cf. a^, = 1/3 
for the standard accretion disk (Shakura & Sunyaev 
1973) and a^, = —0.44 from the composite spectrum of 
SDSS QSOs (Vanden Berk et al. 2001)]. The contribu¬ 
tion estimated by extrapolating the radio emission for 
B1422-I-231 is negligibly small as followed from Figure 2, 
and PG1115-I-080 is radio-quiet. 

For PG1115-I-080, the effects of mid-infrared flux from 
a compact region are also invalidated from the following 
argument. As a working hypothesis, we suppose that in¬ 
trinsic mid-infrared flux prior to lens magnification, de¬ 
noted as fjE, arises both from a surrounding hot dust 
torus, sfjj^, and from a compact region such as an ac¬ 
cretion disk, (1 — s)fjE, where s (0 < s < 1) yields 
the former fraction of the dust torus in intrinsic mid- 
infrared flux. Taking a pair of (Al, A2) images for the 
flux ratio, suppose that a star or a low-mass GDM sub¬ 
halo causes additional magnification for a compact region 
with a factor of {i = Al for Al image or i = A2 
for A2 image), such that a total magnification factor is 
given as (1 -I- where /io,i denotes the magnifi¬ 

cation without a lens perturber. It is also supposed that 
mid-infrared flux from a dust torus remains unchanged. 
Then, the mid-infrared flux ratio, rm, is given as, 

s-I-(1-I-(j2|^)(l - s) 

riR - rsmooth X ■ (2) 

If mid-infrared emission is supplied only from the com¬ 
pact region such as an accretion disk (s = 0), we ob¬ 
tain riR = ropt = ^ smooth (1 + + ^Ai)^ which 
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is inconsistent with the current observational result of 
rjR ~ rsmooth > Topt- Similarly, a more massive sub¬ 
structure with M being large enough to affect mid- 
infrared flux from a dust torus yields rm = Vopt as well 
and is thus unlikely. Instead, if mid-infrared emission 
originates from a dust torus alone (s = 1), we reproduce 
the observational result of rjfi ~ rsmooth- Thus, thermal 
emission of the dust torus is a main contributor to the 
observed mid-infrared flux in PGl 115-1-080, which is con¬ 
sistent with the infrared-bump feature of its SED, and 
the lens perturber that produces the anomalous flux ra¬ 
tios in the optical bands of PG1115-I-080 is either a star 
with sub-solar mass or a low-mass CDM subhalo. 

For B1422-I-231, contribution of a compact region to 
the observed mid-infrared flux is estimated as about 0.2, 
and this implies, as shown below, that a lens perturber 
ought to be more massive than a star with sub-solar 
mass. For instance, we consider the flux ratio A/B, 
whose value in radio or mid-infrared ranging 0.9 —1.0 dif¬ 
fers largely from the smooth lens prediction of 0.76, and 
this anomaly mainly causes the deviation from the cusp- 
caustic relation, (A-|-C)/B= 1. Then, if we set s ~ 0.8 
and if a star affects either A or B image (i.e., = 0 

or = 0), we obtain rm < 0.8 from equation (3) for 
either case as Vopt is smaller than 1 (Table 4). This is 
inconsistent with rjn — 0.94 at 3 cr level. Also, in this 
microlensing hypothesis, the observed rm can be repro¬ 
duced only when s = 0.4 — 0.5 (50 — 60 % contribution 
from a compact region to the mid-infrared flux), which is 
very unlikely. For s = 1, microlensing hypothesis is ruled 
out as well unless otherwise we expect rm = rsmooth- 
Curiously, observations in short wavebands [F480LP and 
F342W (Impey et al. 1996) and also Gunn r and g (Yee 
& Ellingson 1994)] yield the similarity of A/B to the 
smooth lens prediction, whereas the CASTLES data in 
near-infrared band appear to be consistent with the radio 
and mid-infrared results (Table 4). This might be caused 
by differential dust extinction and/or spatial variation of 
QSO absorption medium in the foreground lens but the 
firm conclusion is yet unavailable. 

4.2.2. Effects of Nuclear Starburst and Host Galaxy 

If the mid-infrared fluxes originate from extended re¬ 
gions, such as a region of starburst and host galaxy itself, 
then the observed images would be largely deformed, pos¬ 
sibly in an arc-like form, by foreground lens galaxies; the 
deformation of images from a circular shape is at most 
by a factor of leasing magnification in the tangential di¬ 
rection. This conjecture yields upper limits on the sizes 
of mid-infrared sources. The observed leased images of 
both targets were found to be a circular shape within 
the range of their FWHM, 0."38. Then, taking into ac¬ 
count the largest magnification factors (13.5 for image 
Al of PGII15-f080, 9.9 for image B of B1422-k23I) ob¬ 
tained in the lens model of Chiba (2002), the source size 
Rs for the observed mid-infrared flux must be less than 
~ 240 pc for PGI115-f080 and ~ 280 pc for BI422-t231. 

Effects of starburst activity associated with QSOs can 
be assessed by the excess of far-infrared dust radiation 
as obtained from millimeter and/or submillimeter obser¬ 
vations (Carilli et al. 2001; Omont et al. 2001). For 
PGl 115-1-080, such an excess is unidentified (Barvainis 
& Ivison 2002) as shown in Figure 2. Then, we estimate 


the possible contribution of nuclear star formation activ¬ 
ity to the observed mid-infrared flux using the polycyclic 
aromatic hydrocarbon (PAH) emission feature at 3.3 /rm. 
It is known as an indicator of star formation activity, and 
is located near the rest-frame wavelength of 4.3 /rm for 
the observed 11.7 /im. Based on the reported equiva¬ 
lent widths of the PAH emission for nuclear regions of 
Seyfert 1 (< 1 — 10 nm: Imanishi & Wada 2004) and 
starburst galaxies (~ 100 nm: Imanish & Dudley 2000), 
it is found that the contribution of star formation activ¬ 
ity to the currently observed flux is only 0.01 — 0.1 or 
less, thereby being unimportant for our results. On the 
other hand, for B1422-I-231, no millimeter or submillime¬ 
ter observations have been conducted, so the possibility 
of nuclear starburst in this system is yet to be settled. 
If additional mid-infrared flux from nuclear starburst is 
included in the observed flux, then the anomaly of its 
flux ratio as reported here would be caused by a more 
massive substructure than that based on the dust torus 
size, because the nuclear starburst region would be ex¬ 
tended on a scale of ~ 100 pc or more. In that case, the 
possibility of a GDM subhalo to induce the anomalous 
flux ratio is made more robust. 

Effects of host galaxies appear to be minor, as judged 
from the smallness of the size bounds described above. 
Using the results of the HST imaging observation for host 
galaxies of QSOs at z ~ 0.4 (Floyd et al. 2004) and the 
color information for low-redshift QSOs’ hosts (Jahnke, 
Kuhlbrodt, & Wisotzki 2004; Hunt, Giovanardi, & Helou 
2002), it is found that the contribution of host galaxies to 
the observed mid-infrared fluxes are only less than 0.1 %. 

4.3. Substructure Sensing in PG1115+080 

As argued in §4.1, the observed mid-infrared flux ratio 
A2/A1~ 1 for PG1115-I-080 implies that its anomalous 
flux ratio in optical and near-infrared bands (Table 3) 
may be induced by either a star with sub-solar mass or 
a low-mass CDM subhalo with < 2.2 x IQ'^Mq if 

modeled as an SIS; our current observations do not dis¬ 
tinguish between both possibilities. A possible method to 
settle this issue is to consider the rapid time variability of 
these images caused by microlensing if a star is involved 
in flux anomaly. The time variability of lensed QSO 
emissions can be induced by relative motion between the 
source and lens (e.g., Schneider, Ehlers, & Falco 1992). 
For the purpose of simplification, we assume that this 
relative motion is dominated by transverse velocity (r^j,) 
of a substructure lens, i.e., the sum of its peculiar veloc¬ 
ity inside a galaxy and the motion of the galaxy. Tonry 
(1998) showed that the leasing galaxy of PGI1I5-I-080 
has the internal velocity dispersion of 281 ± 25 km s“^ 
and that the galaxy belongs to a group with a disper¬ 
sion of 326 km s“^, so for simplicity we adopt v_l ~ 600 
km s“^ for the following estimates. 

A useful measure for the time variability is the time re¬ 
quired for a star to cross its Einstein radius with velocity 
v±. For PG1115-I-080, we obtain 

i /I , -.DosOe Dol 

tE = {l + ZL) - - - 

v± Dos 

where Dql and Dos are the angular diameter distances 
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to the lens and source, respectively, and the second equa¬ 
tion is derived for the point-mass lens model. Another 
characteristic timescale is the time required for a caustic 
to cross over a source size Rs- 


tc = 2(1 -I- zl) 


13.8 


Rs Dql 
V± Dos 
Rs 

1-3 


n_L 


6 X 10-3 pc/ V600km s 




yr . (4) 


where we use the estimated source size of an optical con¬ 
tinuum emission region of PG1115-I-080 using the F160W 
observation. Note that the size of a UV continuum emis¬ 
sion region is expected to be smaller by some factor as 
argued in Appendix, so the timescale can be shorter when 
observed in shorter wavebands. 

Although more realistic estimates must consider a com¬ 
plicated caustic network involved in a star field rather 
than a single star, the fact that the optical flux ratio for 
A1 and A2 images of PGl 115-1-080 has remained basi¬ 
cally unchanged over the past decade (Table 3) seems to 
be consistent with the implications from the above rough 
estimates, i.e, the time variability induced by microlens- 
ing is yet difficult to assess because of its rather long du¬ 
ration. However, somewhat shorter and low-amplitude 
variation in these images would be yet possible, espe¬ 
cially in shorter observing wavebands, owing to the ef¬ 
fects of a caustic network in the microlensing hypothe¬ 
sis. In contrast, the time variability is extremely long 
in the millilensing hypothesis and even for less massive 
CDM subhalos only a single event of flux variation may 
be detected because of their much smaller number than 
stars. Thus, tighter limits for substructure in the lens of 
PG1115-I-080 may be set by a long-term monitoring of 
Al and A2 images or that of detailed spectral features of 
their emission lines (Narasimha & Srianand 2000; Mous- 
takas & Metcalf 2003). 


5. CONCLUSION 

We have presented mid-infrared imaging at 11.7 
^m for the quadruple lens systems, PG1115-I-080 and 
B1422-I-231, which are characterized by the anomalous 
flux ratios among their lensed images in other wavebands. 


The lensed images of these QSOs have successfully been 
separated from each other owing to the diffraction- 
limited resolution of the COMIGS/Subaru images with 
FWHM = 0."38. Our calibration of the mid-infrared 
fluxes and flux ratios among different lensed images have 
revealed that (1) the observed mid-infrared fluxes provide 
clear infrared bumps in the SEDs of the target QSOs, in¬ 
dicating that the observed mid-infrared fluxes originate 
from a hot dust torus around a QSO nucleus, (2) the 
mid-infrared flux ratio A2/A1 of PG1115-I-080 is virtu¬ 
ally consistent with smooth lens models but inconsistent 
with the optical flux ratio, (3) the mid-infrared flux ratios 
among (A, B, C) of B1422-I-231 are in good agreement 
with the radio flux ratios. Based on the possible dimen¬ 
sions of dust tori around these QSOs estimated from the 
relation of dust reverberation mapping, these observa¬ 
tional results imply that the origin of anomalous optical 
or radio flux ratios for these QSOs may be caused by 
either of a star with sub-solar mass or a low-mass CDM 
subhalo with an SIS model of Mfjp < 2.2 x IO^Mq inside 
radius of 100 pc for PGl 115-1-080 and a CDM subhalo 
with M^qq > 7.4 X IO^Mq for B1422-I-231, respectively. 

Assessing the presence of many dark subhalos in a 
galaxy is a difficult task because of their dark nature 
but is important for understanding their crucial roles in 
the formation processes of visible parts of the Universe 
as well as for testing the currently standard CDM sce¬ 
nario itself. Lensed QSOs with anomalous flux ratios 
have been, so far, possible sole sites for the investiga¬ 
tion of substructure in their foreground lenses. Further 
mid-infrared observations of other such lenses, especially 
radio-quiet QSOs, would provide useful insight into the 
nature of substructure leasing. 
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APPENDIX 


THE SOURCE SIZE OF AN OPTICAL CONTINUUM EMITTING REGION 

Since the UV and optical continuum emission of AGNs are considered to originate from the accretion disk around 
the central black hole, we estimate the source size by a simple model of the accretion disk. Assuming a temperature 
profile as T (X i?“3/4 jg j-adius of an annulus of the disk), which is expected from the standard accretion disk model 
or the irradiation disk model illuminated by a central source, the disk spectrum is given by 


u 


11.2 Jy 


R 


light days 



(Al) 


which is derived based on the equation (2) of Collier et al. (1999). D is the distance to the AGN and A is a parameter 
determined by a disk model, estimated as A < 4 by Collier et al. (1999). i is the disk inclination angle; cosi < 1 
is expected according to the nnified scheme (e.g. Antonucci 1993) because both target QSOs are type 1 AGNs and 
face-on view is expected. Thus ignoring parameters of A and i, the relation between the wavelength-dependent source 
radius and the luminosity is derived as 


/ XLx{0.51fj,m) \ / A \ 

\ 10^3 ergs s“3 / v^m/ 


R{X) = 2.7 X 10-2 pc 


4/3 


(A2) 
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using fi, oc as calculated from a temperature profile of T oc i? Accroding to this equation, the source size 
was estimated from the optical luminosity and the rest-frame wavelength. 
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